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ie Introduction 

Shock waves have been studied in the past in order to provide insight to the 
energy released during an external shock to a system and the resultant physical and 
chemical reactions which occur. A detonation may be defined as a supersonic wave 
which propagates throughout a crystal lattice structure, and whose energy is released 
during an exothermic chemical reaction.! It is convenient to classify a detonaticn into 
two categories, the first category being defined as the shock front which is 
characterized by a sharp rise in energy and agitation over time occuring on the order of 
picoseconds. As the shock front passes through a lattice, an almost instantaneous 
change in atom position and momentum may be observed. The second category is the 
shock wave which is usually characterized by the slowly decaying nature of the wave 
until thermal equilibrium has been established.” The shock wave transfers energy 
through the lattice in a much slower fashion. Time delavs are experienced with respect 
to atom displacement as well as energy transfer. Shock waves propagate through most 
materials on the order of 2500 meters per second ( 5.49 x 10° a,/Ps ).1 There has 
been considerable attention in the past devoted to the examination of shock waves and 
whether they may be classified into two seperate categories ( first and second sound ). 
This idea has been disputed in the literature for several years.3-5 

Shock waves produce a variety of anomolous optical and electrical properties as 
well as nonthermal chemistry, and they are capable of producing bond dissociation, 
phase changes, and a redistribution of vibrational, electronic and rotational energy 
within a lattice. Shock waves are also unique in that energy may be deposited to 


vibrational modes of 2 molecule under nonequilibrium conditions. 





Early studies concerning simulations of shock waves propagating within crystal 
lattices were restricted to classical Gynamics with finite dimensions, consisting of non- 
oscillatory atoms prior to the introduction of a shock wave.32+ These studies were 
primarily aimed at understanding the shock wave, and its velocity and kinetic energy 
effects upon a crystal lattice. Later studies>-8 devoted considerable attention to the 
characteristics of the shock wave and the various phases of velocity transformations. 
Classical studies have been completed on systems which were initially at a nonzero 
temperature (zero temperature being defined as no oscilatory motion),?-!9 and on 
models which incorporated deformations within the crystal lattice.1!-14 These studies 
were concerned with affects of free radicals, mass defects which were randomly 
dispersed, and mass defects of heavy and lignt impurities on shock wave propagation 
throughout the lattice. 

There have been few quantum mechanical studies concemming shock waves. 
Dancz and Rice derived expressions for the quantum mechanical equations of motion 
for coupled anharmonic oscillators.!> There have been other studies which were also 
devoted to deriving equations of motion for quantum systems by the introduction of 
raising and lowering operators for interacting Morse potentials.!© Other studies 
examined the probability of occupation of excited states and the state to state energy 
flow in which energy became trapped in excited states.!/ 

Experimental work using picosecond lasers has also been done in order to 
understand shock wave propagation through lattices. In this experiment, a 


gl2 


picosecond laser ( capable of 1 Watts/em* ) imposed a shock upon a sample of 


water which was analyzed using Raman Spectroscopy. This was the first experiment 





in which energy transfer from a macroscopic shock to intermolecular states has been 
observed. . 

Most of the previous studies have been devoted extensively to the classical study 
of shock waves, and the subsequent interface with a lattice model using classical 
mechanics. Much work has been devoted to examining the shock wave and the 
resultant energy transfer associated with it. Recent studies [Wyatt and Marston] have 
been devoted to describing a shock wave effect quantum mechanically in a one 
dimensional crystal lattice. This work was primarily devoted towards a high energy 
shock wave and its effect upon a completely harmonic system. A classical linear model 
of twenty diatoms was formulated and perturbed via a ballistic particle of equivalent 
mass. The inner two diatom pairs were chosen for the classical study in order to 
restrict center of mass movement within the lattice and to represent inner lattice 
substituents. These inner atoms were perturbed by direct interactions of their neighbor 
atoms (entering and absorbing driving atoms). Wyatt and Marston used the action of 
these atoms to establish a time-dependent driving potential for two diatoms. The 
entering atom (q<) provided the initial interaction from the shock to the cluster; 
whereas, the absorbing atom (q>) acted as a resevoir for the energy to be released from 
the cluster. The dynamical data from the classical results of qe and q> was used 
directly to formulate an interaction potential for a quantum mechanical model of two 
diatoms; thus, the classical and quantum model experienced similar potentials. The 
quantum model was represented by four normal modes, upon perturbation yielded 
quantum dynamical data for comparison to the classical system. Wyatt and Marston 
showed that a high energy shock imposed upon a harmonic quantum system 


converges, as expected, to the classical results. 





This work is devoted to the replacement of harmonic oscillators with cubic 
anharmonic oscillators and in making comparisons to classical dynamics. A similar 
anharmonic classical model is established in order to provide for a driving potential in 
our quantum svstem. In the quantum system the cubic anharmonicities provide vanous 
mode couplings; thus, our Hamiltonian may be not separated exclusively into the four 
modes as in Marston's; however, it will be represented by all four modes. This fact 
will restrict our quantum basis size. 

The remainder of this thesis is organized in the following fashion. Section II 
develops the model for both classical and quantum systems. The classical equations of 
motion are derived showing how we obtained dynamical data for comparison and for 
the quantum potenual. The quantum coupled anharmonic Hamiltonian is derived using 
the normal modes representation and the tme dependent svstem is developed by 
approximately solving the Schrédinger ume dependent equation. Quantum formulanon 
for posinion expectation values, normal mode energies, and bond energies 1s presented 
in this section. 

Section II presents results of quantum harmonic and anharmonic calculanons 
of several basis sizes and compares these results to classical dynamical data. Section 
IV draws conclusions based upon thse results and also provides commentary based 


upon different formulations of the problem. 
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Ie Method 

Model 

A one-dimensional twenty-diatom classical model was established in order to 
study the effects of the propagation of a shock wave through a crystal lattice. All of 
the atoms in the chain are of equal mass (mass of fluorine atom); the intermolecular 
distances represent the intermolecular separations of diatomic fluorine. Each diatom 
was connected with a force constant Ks and equilibrium separation po. The 
equilibnum separation between diatoms was Rog with a connecting force constant Ke 
(see Figure 1 and Table 1). A ballistic particle of equal mass will impose a shock wave 
upon the chain of atoms. In order to compare a classical model with a quantum model, 
the central pair of diatoms in the classical one dimensional chain will be isolated. This 
pair of diatoms is far enough from the ballistic particle to reduce any sporadic results 
from the strong potential and a smaller system of two diatoms was chosen for 
comparison to an equivalent quantum model based on the fact that a twenty diatom 
cluster is numerically intractable to handle quantum mechanically. Classically, the 
central diatom pair will be perturbed via interaction with the nearest atoms in the chain. 
The entering atom (q<) provides the initial interaction from the shock to the cluster; 
whereas, the absorbing atom (q;) acts as a reservoir for the energy to be released from 
the cluster. The dynamical] data from the classical results of q-(t) and q(t) will be used 
directly to formulate an interaction potential for a quantum mechanical model of two 
diatoms; thus, the classical and quantum models will be experiencing similar 
potentials, thus allowing for direct comparisons. 

Classically, these atoms may be thought of as forty harmonic or anharmonic 
oscillators. At time t=0, atom positions were arbitrarily set at their equilibrium values, 


5 





and atom number one in the chain was originally chosen to have zero momentum. 
Atom momenta aN easily obtainable through the iellowtin relationships: !8 

PEEP TKS (G2- 4 -p.)+1G2-41-P,) = Es (1) 

Pi = momentum of atom number one = 0 

Oo = equilibrium separation (Table 1) 

(q2 - qi) = Po; since at equilibrium 

Ey = Total energy (sum of potential, and kinetic terms) 

tu == reduced mass (Table 1) 

y = degree of anharmonicity (Table 1) 

Ks = Force constant between atoms (Table 1) 
The total energy (ET) between any two atoms was calculated from the ground state 
energy of a diatom pair (ET = ae where @ =V K/u1). Equation ( 1 ) readily reduces to 
Oe a 
on ~ “T , since the potential energy term and the anharmonic term consistently drop 
out of the initial equations because all interatomic distances are at their equilibrium 


separation. From this equation, the momentum of atom number two can be calculated, 


thus allowing all other momenta to be calculated for their initial values via: 


2 
(Piet fi: Pr) = Ee 
2u 


In order to initiate a shock wave into the forty diatom chain, a ballistic particle 
(impact atom) with equal mass of a flourine atom imparts a potential upon the system. 
The interaction potential of the impact atom and with atom number one 1s given by the 


following equation: 


Mee. Os - 4, 
Vac 
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a = 1.0(Bohr!) 

q) = position of chain atom number one 

do = position of impact atom 
The quantity A was chosen as twice the ground state energy value of a diatom pair, and 
a was chosen as 1.0 Bohr -! to optimize the trajectory of the impact atom in order to 
produce a movement of the atoms at a speed greater than 2500 m/sec (speed of a shock 
wave). The entire classical potential of the chain may now be wnitten as a combination 
of the interaction potential from the ballistic particle, the sum of harmonic potential 


interaction of the particles and an anharmonic term: 


40 40 
-a(q,- 1 2 l 2 
V=Ae = eee oe > Ke n+i- In- Po) as ys ae peer 1s, ) 
Nees oa n=2,4,6..._ 
40 3 40 3 
+ » ieee eee Ds) 17 Dy COL aay = el ss 
Wade ose n=2,4,6... 


where 
Gn = position of atom n 
Oo = intramolecular equilibrium distance (Table 1) 


Ro = intermolecular equilibrium distance (Table 1) 


As the ballistic particle moves toward atom one with twice the initial momentum 
of atom number two, the repulsive potential becomes more significant, thus creating 
transverse motion within the chain. These atoms move according to Hamilton's 


equations of motion: 


nS eae eo | eel, 
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Upon numerical integration, these coupled differential equations yield atom momenta 
and positions as a function of ttme. The Adams-Moulton integrator was used in this 
study with the following parameters: 

Order = 6 

DT = 1 atomic time unit (2.5x10-5 ps) 

EPS = 10-3, maximum relative degree of iteration 

N = 60,000 = numbers of time step iterations 
Results produce classical trajectory values for the one dimensional system. As stated 
previously, the two diatoms located at the center of the cluster are of main interest for 
this classical problem and for comparison to the quantum mechanical problem. 

In the quantum mechanical system, the four atoms are initially all located at their 
average equilibrium value corresponding to their counterparts in the classical system. 
These four atoms are connected to driving atoms (q< entering atom; qs absorbing atom) 
as formulated in the classical problem, and the motion of these driving atoms in the 
classical system is converted into the driving potential for the quantum system. This 
potential creates a perturbation within the quanutm model causing atom displacement. 
Since the one dimensional quantum mechanical system may be thought of as four 
harmonic or anharmonic coupled oscillators, motion in the system 1s described by four 
normal modes. The positions in the system may be expressed as a sum of initial 
Cartesian position plus a time dependent displacement 6: 

qi=4i +8 (2a) 

q2 = 42+ 89 (2b) 


re) 
2 — 3 03 (2c) 





O 
q4= 44+ 04 (2d) 
From these relationships, calculations of transverse atom movement are performed. 


Formulation of the Hamiltonian for the Cluster 
For a set of atoms which are interconnected via sets of anharmonic springs, 
there are many types of internal motions, vibrations, which are dependent upon the 
initial displacements of the atoms. The Hamiltonian for the anharmonic system may be 
expressed as the sum of a harmonic part (H°), anharmonic part (Va), and tme 
dependent part (V(t). 
Hy =H° + V, + V(t) 


Formulation of Harmonic Hamiltonian 
The harmonic Hamiltonian for this one dimensional two diatom case 1S 


expressed as the sum of the potential and kinetic operators!? 


H° “SP sa qe(t)-R.)° + qs (0 -q4-Ro)‘] 


Ky 2 2 Kee 5 
+ [(d2-41- Po) + 4-93- Po) 1+5- (43 -92- Ro) 
where 
p; =m, q, = momentum of atom i 


mj = mass of atom 1 


qj = position of atom 1 


Oo = equilibrium intramolecular distance 
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Ry = equilibrium intermolecular distance 
Kz = Force constant of intermolecular bond 


> = Force constant of intramolecular bond 


As stated previously, the positions of the four atoms as a function of time may 


be described as the sum of an initial position and a displacement position (6). The 
positions of the impact and absorbing atoms (q< and qs) may also be expressed as an 


initial position plus a displacement: 


qc(t) = qe + S(t) (3a) 
qs(t) = 45 + S(t) (3b) 


The harmonic Hamiltonian is time independent, thus d<(t) and d5(t) are equal to zero 
until a perturbation is introduced. Utilizing the previous relations of initial positions of 


atoms, it is clear that: 


Oo oO 
Go- Gi = p- 
0 0 

iA Ciat== 1) 
q3-q,=R, 
q.-q,=R, 
gi GS Se 


This enables simplification of the harmonic Hamiltonian: 





4 
wD SS ( oe ss as as ycls =? (853,)" (4) 


Formulation of Harmonic Time Independent Hamiltonian Matrix 


If an atom has an original position g° and throughout time moves to a different 


position q(t) the resulting vector describes the displacement (6) from its equilibrium 


value . This may be represented in mass scaled displacement coordinates (1;) 


ni = (mi) 1/25) 


In matrix notation, 1) and 8 are defined as a row vector, and iv 


fil/2 is a square matrix 
with the root values of atom masses on the diagonal 


The kinetic energy may now be represented as 


or in matrix notation 


where 7] ‘is the transpose of the row vector 7}. 


I, 





The potential energy of the system may also be represented in the same manner 


We expand V in a Taylor series in the mass scaled displacements: 





a nye yet P Enns 


j=l 


where 


V = potential energy 
Wie v (0=0) =0 


yn rads =Q0 (Expanded about 6=0, the minimun of V ) 
| 


Thus, the harmonic potential energy may be expressed as 


2 
! dV os 
Vaz) 20; ) = Kn 
- - a | 2 
ej onan, 
where 7} 1s the mass scaled coordinate row vector, Tj! is its transpose (a column 


vector), and K = force constant matrix which may be represented as follows: 


K 32 
S| eK “Ke 0 
POs K 2 Ke-Ko 

Opn On Ke K 


Ke is the weak intramolecular force constant 





Ky 1s the strong intermolecular force constant 
K =K, + Ke 
In order to analyze the system, an orthogonal transformation is required such 


that At'=A-!. Define € =A and C= A'yt where C is the normal coordinate row 


vector. The kinetic and potential matrices may now be rewritten with this substitution. 


T=LC(A N(ANG) = EC 
ce §SMINZ2A d = CAM 1/2 
= (5 )C(A KA)C ‘ = = (5)NA(A KAA] ' 


The AUXA matrix may be redefined as the diagonal normal mode frequency square 


matrix (£22). The potential energy is then 


l 
Wes 56O2Ct 
From these relations the normal mode frequencies have been established (Table 
2). The displacements of the atoms from their respective equilibrium positions may be 


represented as linear combinations of the normal modes 


§ = CAIM-1/2 


where a representation of the A matrix follows: 


epee. 0.9 —-0.5 
eee -Os5 05°05 
eV) -0.5° -0.5 
som. 0.5. 0.5 
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During the analysis of the two diatom case the values of the A matrix were all 
approximately +0.5 or -0.5. This enabled the € operators to be presented as simplified 
combinations of the displacements of atoms from their respective equilibrium positions. 


(pec hig for description of modes. ) 


Cy, =—0.5Vm(5; + 52 + 53 + 54) (5a) 
Co = +0.5Vm(-5; - 52 + 53 + 54) (5b) 
C3 = +0.5Vm(8) - 82 - 83 + 84) (5c) 
C4 = 4+0.5Vm(-5) + 5 - 53 + 84) (5d) 


Formulation of Anharmonic Time Independent Hamiltonian 
If one considers an arbitrary potential V(q) between atoms n and n+1 which has 
a minimum at q=qpo, the potential function V(q) may be expanded in a Taylor series 


about qg:29 


V(q)=a@ +B(qn+1 = hers do)? + ¥(dn+1 - Gn- do)? 


Gamma (¥) represents the physical size of the anharmonicity (deviation from harmonic 
potential well ) and is negative. If one were to make the assumption that the (Gn+1 - dn 
- dg)? term is sufficiently small compared to the first terms of the expansion, a 
harmonic oscillator would define the system. The (qn+1 - Gn - do)? interaction of the 
potential function cannot be ignored in real systems; thus, this term has been added to 
the total Hamiltonian for the model. In order to generate actual physical models, these 


cubic anharmonicities between the two diatoms and interconnecting sets of atoms were 





i 


introduced. The anharmonicities for the two molecules in the cluster are (54 - 63)3 and 
(52 - 1)>. We are assuming that anharmonicities exist only between diatom pairs, and 
that only the cubic anharmonicity is significant. 

From equations (5), the following normal mode operator equations were 
derived: 

“(54 = 63)3 = ym (C, cee F 

y(S2-8)3=7M (Cyt 03) 

Anharmonicities within the intramolecular bonds couple only modes three and 

four. If the intermolecular bond 1s anharmonic, all four modes would couple. 


Expansion of these terms resulted in the following anharmonic potential operator: 


Sf 


ae 2 
Va=ym (20, + 663 64) (6) 
Formulation of Time Dependent Hamiltonian (V(t)) 

In the time independent Hamiltonian, the impact and absorbing atoms were 
stationary; thus, the d<(t) and ds(t) were equal to zero. Now movement of these atoms 
must take place in order to propagate a shock through the cluster. The terms two of the 


harmonic Hamiltonian containing these groups were as follows: 


Ke Oo Oo Z Kz 2 
z (41-454 91 - o_(t)-R,) Sa = ON) 


Ie fe) fe) Ke 2 
SS (a3- a4- 84438, - Ry) = 556, - 8) 


Expansion of these term yields: 





1 6 


IS : 
=z (81? - 28) 8. + 8.2(0) + 85(0)? - 28,844 842] @ 


The time independent terms, (517, 542 ), are formulated in the time independent 
Harmonic Hamiltonian Eq. ( 4), thus equation (7 ) may be seperated into time 
dependent and time independent terms. The time dependent terms will produce our 
driving potential for the atoms in the cluster. Expansion of 5; and 64 in terms of 


normal mode operators (eqns 5) C's yields: 
Se ~12 =12 
pent OC = 3 Cay tom «(Cy > C203 — C4) + &.(t)? + 82(t)] 


Movement of the impact and absorbing atoms therefore generates a change in position 


which results in a potential being developed. 


Definition of Raising and Lowering Operators 
The operators which correspond to momentum and coordinates must satisfy the 
commutation relation. In matrix mechanics, the one dimensional harmonic oscillator 


Hamiltonian is as follows: 





La 


In order to procced with normal mode analysis, the following raising (at) and lowering 
operators (a) are introduced.!? . 

a= apr + BC 

at= a*pe + Bt 
c*B is chosen as a purely imaginary value such that ap = —o 


77 


C= 


+ Z 22 1 
a a= [pp + Q C ies 
ee 
— 
— = (+ipc + QD) (Sa) 
lige: od : 
a= ae (—Ipc + OC) (5b) 
So the Hamiltonian of the system may be described as 
hOQ 
H =a "a 7. me 


The Hamiltonian may now be applied to the Schrédinger equation in order to evaluate 


the eigenstates. 


wee E ES 
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@'aa l= Ele 
‘ale =(E-**) | E> 


such that 
AQ 
<Ela'‘a se (E-——)<E | E> 


28 ES Si 


AQ 
Sela ale>=(b——_) 


This will provide a harmonic basis for the problem. The states may be labeled as 
numbers with corresponding energies evaluated from the numbered states. 
The raising and lowering operators may now be defined through their action 


upon the eigenstates. 


aln>=~VnhQ | n-1> (Ya) 


a! | n> = V (nt+)hQ | n+l> (9a) 


In the model problem, the normal mode operators were defined from equations ( 8 ): 


C=(Aza) ata’) (10a) 
pe=(/ys) (a-a) oy 
Using the above relations for a and at the cubic anharmonic operators may now be 


evaluated 


he nearest) nee 15 


CiIn> = ( 








eS, 








2 A 
ee ret a oeeesnd |n, > 4 (n:4 1) ln, > 
m;Q. - 
mc ee 2) net 25) 
3 3/2 1/2 1/2 2 2 
C. Ine =) on ecreea( ee) — \meaSecr(n) (rei ce(ay. (lee 1) 


an | 
iS 6 
++ (n;) ine (n; + 1) es In; - I>+ (yaoi) ae (n, +1) (n; +2)) In. + |> 


ee (n.2). (ne? 3) > Inet 35] 


The original basis set chosen for the problem consisted of the product of four 
normal modes ( | n\> | n2> | n3> | n4>). The normal mode operator acts only upon its 
corresponding basis state (i.e. C43 operated only upon | n4>). The normal mode 
operators were used to formulate the anharmonic and time dependent Hamiltonian 
matrix elements. In the total Hamiltonian matrix, the harmonic terms correspond to the 
diagonal elements and the anharmonic terms correspond to the symmetric off-diagonal 


elements. 


<n; IC. In> =( vi (Ga ae aes 





ihe Ni Nis) 


Mj +6; 


The anharmonic Hamiltonian matrix was generated and diagonalized to obtain the 
appropriate eigenvalues and eigenstates ( Table 3 ). 
Topol ele 


— dod 


where Z = eigenvector matrix 





Z' = transpose of eigenvector matrix 


E = diagonal eigenvalue matrix 


Propagation of Time Dependent Hamiltonian 

The following passages formulate the tme dependent perturbation introduced 
into the one-dimensional system.2! The basic solution to the time dependent 
Schrédinger equation is as follows: 


es 
Hy = Gh) 


where H is a function of position and time in a Cartesian coordinate system. This 1s 
extremely difficult to solve; therefore, the Hamiltonian is divided into a time- 
independent and time dependent part: 

Bie i, + Y(t) 

Hy] = time independent tenn of Hamiltonian 

V(t) = time dependent term of Hamiltonian 

PY = (tye KEM 
The spatial equation to solve is 

Hry = Ew 

E = separation constant 
It has well behaved solutions when E is equal to En 

AryWn°? = Enyn? 
Assume the ‘¥(q,t) can be written as a linear combination of orthogonal functions: 


Sa —= Sean qv) 


O oO : 
_ -E ¢/h 
Y =e 
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oO 
The quantity ‘V’ , is described by the eigenfunctions of a harmonic oscillator and the 


quantities ap(t) are defined as the occupation coefficients. Since the wave function is 


normalized, 


> a, (0*a, (0 = 1 


the series expansion for ‘¥(q,t) may now be substituted into the Schrédinger wave 


equation yielding: 


d h = O* oO 
Fant) =-+ >, an i We, VOW, dq 


n= | 


witemeanmae= I, 2. 3, a. 
This yields the rate of change (with respect to time) of the m“ component of the 
basis. Am (rate of amplitude change) depends upon all basis amplitudes. The 


following equations are shown to exemplify the above differential equation. 


h da, o* ° o* o ae o 
See 1) Vs V@w,dt+ta.} yw, VQwjdt+...tay | w, Vi) wydt+... 


A, da oF fe) OG 6) oF 6) 
Saas fv; vowara fv; VO Wadi notag f Wo V(t) Wy dtt+ ... 


Utilizing these coupled first order differential equations, numerical analysis using the 


Adams Moulton integrator provided the coefficients for the various states of the model. 
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Calculation of Quantum Position of Atoms 
In order to calculate positions of atoms, computations of expectation values of 
displacements from equilibrium were performed. The matrix elements of operators 6; 
were evaluated in the normal mode basis, and then transformed into the harmonic 
representation using the eigenvector matrices , as follows: 
Z'D.Z=D. 
where 
Z' = transpose of eigenvector matrix 
“= eigenvector matrix 
[Dj = displacement matrix 
1D'"| = transformed displacement matrix 
Upon transformation, D'j was then summed over each state with its corresponding 
energy and probability: 
PQ) D. P(t) 
where 
[P(t)]i=ci(Qexp[-iEjt/h] 
where P is a column vector based on probablity and eigenvalues, and P! is a row 
vector. This summation yields the displacement of atom1 from its equilibnum 


position. 


Calculation of Quantum Bond Energies 
The total energy of a particular bond is expressed as the summation of the 
bond's potential and kinetic energy. The potential energy for the quantum system 1s 


calculated from the expectation values of position for the corresponding atoms. 





pee: 


Potential energy = (1/2)(Kpond)(9j - 8)? 


From equations (5) the potential energy can be represented in terms of normal modes: 
I 


Z 
PE12 =) (K,) (8,-8))?= ©) (K,) G,- 65) 

Zz 
PE23 = (5) (KQ (83-8)? = ©) (KO (65-69) 


2 
PE34 = (5) (K;) (@4- 83)?= (2) (K) (Ca + 63) 


These operators will act upon our normal mode basis set in order to establish potential 
energy matrices. Upon transformation of the potential energy matrix with the 
eigenvector matrix, summation over all states at a particular time yields the potential 
energy of a bond. 

Z'PE,.Z=PE,, 

P (t) PE P(t) = potential energy of bond y 


As stated previously the kinetic energy 1s calculated by: 


(2 
CSS niA 
T=366' 


where = is a column vector and is normal mode momentum. As mentioned with 
respect to position, normal mode momentum ( Pz; ) may be expressed as a summation 


of corresponding atom momentum (pj) 


Prey = —0.5Vm (p; + pa + P3 + Pa) = (IKI V2) (ay-2 1) 
Peo = +0.5Vm (—p; — p2 + p3 +Pa) = ( 1Ai¥2)) (ay-ay) 





pe3=+0.5Vm (py — p2— p3 + py) = (IMiV2)) (a3-a3) 
Pca= +0.5Vm (~p, + po — p3 + py) = (1/G¥2)) (yay) 


The kinetic energy of the 1-] bond is expressed as: 


eye 


i 
j ay 


which is expressed by normal momentum operators as: 


2 
(Pe4- Pcs) 
T)o= passes kinetic energy operator for bond R;> 
ou 
2 
(Pe2- Pea) 
T53= ei eal aa kinetic energy operator for bond R43 
Zu 
(Peat Pes) 
T34= gee ee’ = kinetic energy operator for bond R3, 


2u 


Each kinetic energy operator established a kinetic energy matrix which 1s transformed to 


the harmonic representation via the eigenvector matrix Z. 


t 4 
cE..Z4=RKRE-p 
Te, Ke,2Z iN Boi 


Upon transformation we sum over all states: 


P () KE; PO) 
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Calculation of Classical Bond Energies 

Classical potential and kinetic bond energies are calculated directly from atom 
position. 

Potential energy = ()(Kipond)(Gn+l - dn - Po)? 


2 
Kinetic energy = Cs. (ome =e) 


Calculation of Atomic Period 
The atomic period for this model was based upon the stonger force constant. 
M=(Ks/u yl? 


T=2n/@m=2.1 x 103 atu 


Energy Evaluations of the Two Diatom Model 
In order to evaluate the normal mode energies, the original Hamiltonian matrix 
was separated into four normal modes, a tme-independent anharmonic term, and the 
time-dependent potential term. 
Hp=H,+H,+H,+H,+V,+ Vi) 
Each matrix was initially transformed into the harmonic basis by multiplication with the 
eigenvector and transpose eigenvector matrix (Z = eigenvector matnx). 
ZH.Z=H. 
Upon transformation, each term was summed over all possible states to provide the 


energies of the particular parts. 


PG) H, P(t) 





III. Discussion of Results 

A. Classical Trajectories 

When the ballistic particle imposes a shock, it adds .005884 Hartrees of energy 
to the chain of atoms. The shock wave propagates through the classical chain at 
approximately 6.5 x 10-2 ap/ps, causing displacement of atoms ( Figures 3,4). The 
classical system responds as expected with respect to atom position and bond energies. 
Atom position 1s directly related to the neighbor atoms, and atoms tend to oscillate in 
diatom pairs. The impact atom provides a disturbance within the lattice which is not felt 
by the central atom pairs for approximately eight periods (Figs. 5). The transverse 
displacement of the entering atom (qo) imparts energy to the first atom in the lattice 
which causes the transverse movement of this atom and so forth; thus, resulting in the 
propagation of the disturbance throughout the lattice. Since this shock wave is a 
comparatively low energy disturbance with respect to previous calculations 
[Wyatt,Marston], atoms are not significantly displaced and atom separation 1s fairly 
constant. The atoms connected by a strong bond tend to move as a pair, with most of 
the lattice deformation occuring within the intermolecular bond. 

Bond energies indicate that energy is transferred between kinetic and potential 
energy . As the energy from the shock wave is transferred to the bonds, 
nonuniformities in the oscillator behavior are observed. The strong bonds exhibit a 
tendency to remain relatively constant, while most of the energy from the shock is 


transferred to the intermolecular bond ( Figures 6-11 ). 
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J65. Convergence of the Quantum Model 

The time dependent coefficients for the quantum states aj (t) will help us in 
examining the convergence of the results of our system, the nature of occupation of 
various quantum States, and the conservation of probability. Quantum eigenvalues of a 
sixteen state model are listed in Table 3. In the present model, the classical positions of 
q< and qs are used in calculating a driving potential. The classical system was 
initialized with all atoms in the chain located at their respective eqilibrium values. Initial 
momenta were calculated for these atoms, and Hamilton's equations of motion were 
solved for position and momentum. Upon solving the equations of motion, the atoms 
were no longer located at equilibnum displacements. This has been defined as t = O for 
the system. Since q< and qs were displaced from equilibria, a small driving potential 
for the system was established causing population of excited states in the quantum 
model. The model may now be thought of as being non-zero temperature. ( atoms are 
oscillatory ). 

When a shock wave is introduced to the quantum model, energy is transferred 
into the crystal lattice resulting in excitation and population of higher energy states. 
This corresponding population of states may be represented as aj*aj where aj 1s the 
coefficient of occupation of state 1 ( .04% deviation in the 16 state model, .0087% 
deviation in the 100 state model). Plots of a 16 state basis illustrate a penodic transfer 
of energy from state to state and a tendency to retain most of the energy within the 
lower mode states (Figures 12-27 ). Harmonic and anharmonic 16 state bases exhibit 
similar tendencies. Examination of a larger basis (100 states) shows a larger retention 
of energy in lower mode states. Since there is a ratio of 5:5:2:2 (nj:n2:n3:n4) in the 100 


state system compared to a 2:2:2:2 ratio in our lower (16 state) model, the energy 





distnbutes itself differently; thus, allowing for even less population of higher energy 


* 


levels (Figures 28-33). 


G; Comparison of Position Trajectories 

Comparison of expectation values for the positions of atoms for several bases 
indicates fairly good convergence of the system. (Figures 34-41 ) We have learned 
from [Wyatt and Marston] that quantum harmonic atom positions at high energy 
converges towards the classical results. Although there are deviations from this 
expected result one can see that the larger harmonic basis sets move toward 
convergence to classical results. (Figures 42-45 ) This is also the case with 
convergence of the anharmonic system. 

Anharmonicities produce deviations 1n energy level spacing as compared with 
a purely harmonic system. The deviations in energy spacing produces a variance in 
energy state occupation; thus, there will be variances in position plots of harmonic and 
anharmonic models. Since there are only two states 1n the higher frequency modes, the 
100 state system does not fully describe the cubic anharmonicity making differences 
between the harmonic and anharmonic models small. (Figures 46-49 ) Atom motion in 
the anharmonic model deviates slightly in frequency and amplitude as compared to the 
harmonic model. 

In comparison to classical harmonic trajectories, the quantum harmonic model 
does not converge exactly to the classical harmonic trajectory, as expected [Wyatt, 
Marston]. The initial oscillatory motion of the quantum model prior to a shock is 
deviates slightly coupled with small deviations of average position away from average 


equilibrium values. This result is most likely caused by the low number of energy 
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states in the quantum model. If one were to examine the frequencies of the four modes 
(Table 2), one can see that they are in a ratio of ieee (205-3207). 

The quantum model currently has a 5:5:2:2 basis ratio. Movement in each 
system should be characterized by oscillations of diatom pairs. For instance, normal 
oscillations of atoms one and two should be directly out of phase as observed in the 
classical system; however, careful observation of quantum diatom pairs yields small 
phase dissimilarities which may be caused by the population of various energy states. 
Since we have a low number of internal mode states (modes one and two), the higher 
modes are increasingly populated. If one were to examine the modes (Figure 2), one 
can see that linear combinations of the modes describe the normal oscillator motion of 
the atoms. Since the higher energy modes are more populated, we observe phase 
dissimilarities within diatom pairs. Classical anharmonic results versus quantum 
anharmonic results correspond in a similar fashion. 

The quantum harmonic model of [Wyatt, Marston] was easily separated into 
rour independent modes. This allowed for vast numbers of low frequency mode states 
as compared to the higher energy modes (500:200:100:100) (nj:n2:n3:n4) giving a total 
basis size of 109. The low frequncy modes were now capable of describing center of 
mass and oscillatory motion. Addition of this number of states to the present model 
makes the problem numerically intractible. 

It is also interesting to note the dependence of the present model upon qs. The 
absence of the potential created by qs causes great phase dissimilarities with respect to 
diatom oscillations, and it also inhibits proper atom atom displacements with respect to 
a shock induced lattice. Several cases have been examined ( Figures 50-59 ) in order to 


show these effects. Since the qs potential is no longer present prior to a shock ( ramp 
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shock ) or not present at all, this causes the excited state population to vary 
significantly; thus, atom energy and displacement peeotie less uniform. The potential 
caused by qs has a significant affect upon lattice dynamics. 

Another interesting point to note is the dependence of this model on the size of 
the anharmonicity. The degree of anharmonicity was increased 10 fold as a trial 
experiment. Small deviations in position are observed as well as small frequency shifts 
in oscillations Fig. ( 60-61 ). In the present model, we are limited to a ground state and 
first excited state in modes three and four. This limitation in basis size effectively 
reduced the survival terms of the anharmonic operators; thus, allowing for small 


deviations in frequency and amplitude of atom positions. 


Be Quantum Mode Energy Comparisons 

A slightly larger difference between our quanutm harmonic and anharmonic 
models is observed in mode energies. (Figures 62-65 ). Energy level spacing in the 
harmonic model is different as compared with the anharmonic model in which energy 
level spacing decreases towards higher levels. This spacing difference exhibits 
variances in mode energies prior to and after shock wave introduction. If both models 
were initially stationary states, one would observe that mode energies would start at 0.5 
quanta. As stated earlier, t =O corresponds to a time in which atom positions are 
displaced from their equilibrium values causing a driving potental to be created by qe 
and qs. This deviation in g<- and qs is fed into our quantum system providing for 
population of some excited states. Since the excited states became populated, initial 
values of mode energies increased above 0.5 quanta. Initially, the harmonic model 


distributes energy to its excited states in a different fashion as compared with the 





anharmonic model. The anharmonic model slowly distributes its energy to higher 
States then reaches a steady value prior to the shock ere. Energy transitions are 
much smoother in the anharmonic model because energy levels have been lowered. 

As shock energy is introduced to both systems, the anharmonic model reacts 
more quickly and distributes energy to higher states. This is evidenced by the fact that 
the higher modes reach a steady state at slightly higher values as compared to the 
harmonic model, and the lower mode energies in the anharmonic case achieve smaller 
energies as compared to the lower modes. In each case, the upper modes pick up 
approximately 0.2 quanta of energy as compared with approximately 0.6 quanta for the 
lower modes. A ratio of this (0.6:0.2) corresponds quite well with a ratio of states 
(2) 

Addition of states to the lower frequency modes in these models will definitely 


improve the energy transfer prior to and during the shock wave. Energy will be 


distribute differently among the lower states allowing for the system to stabilize. 


E. Comparison of Bond Energies 

The classical models react as expected with respect to bond energies. Diatom 
pairs are oscillating in a uniform fashion prior to a shock wave resulting in kinetic and 
potential energy to be directly out of phase and total energy to remain constant. As the 
system is perturbed, oscillations and dispacements vary; thus, variances in potential, 
kinetic and total energy may be observed. The strong bonds have an inherent tendency 


to resist change with respect to a low energy shock. ( Figures 6-11 ) 
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Prior to a shock in the quantum models, energy within the bonds is dependent 
upon the coupling of various modes. Since the system A non-stationary, population of 
various energy levels vary with time. This causes individual bond energies to fluctuate 
transferring energy between bonds. As the shock wave is introduced to the system, 
bonds exchange energy slightly faster and the overall energy within the bonds is higher 
( Figures 70-79 ). There are few differences between the harmonic and anharmonic 
quantum models; however total bond energy is slightly lower in the anharmonic 
system. Careful examination of the potential versus kinetic energy plots of bond 
energy in both quantum systems shows that potential and kinetic energies are in phase. 
Part of this is due to the sharing of energy between bonds; however, some is caused by 
an small number of low frequency mode states (modes one and two ). Excessive 


energy is placed into the higher mode states; thus causing an imbalance in the system. 


This was not as evident in positions as it is in bond energies. 








ae, 


IV. Conclusions 


Examination of a crystal lattice using a coupled mode basis yielded a great deal 


of information concerning both harmonic and anharmonic models. The following 


conclusions may be drawn from these results: 


(1) 


(2) 


(3) 


(4) 


(5) 
(6) 


(7) 


The deviations in energy spacing from the anharmonicities cause 
transitions to occur more easily as compared with the harmonic mode! 
Addition of states to the quantum model allowed for convergence toward 
classical results. 

The anharmonicities produce slight deviations in frequency of oscillation 
and displacement because of energy spacing differences. 

The quantum model is not only dependent on the potential created by qe: 
the potennal created by qs is necessary to allow the atoms to act as a 
cluster. Any deviations in the potential created by qs causes the model to 
change significantly. 

Energy spacing differences allowed for deviations in mode energies. 
Bond energies within the classical models correspond directly to atom 
positions. There are instances in time in which potential energy is at a 
maximum (corresponding to maximum deviation of atom position from 
equilibria) and minimum kinetic energy. Total energy is constant prior to 
a shock. Most of the energy is transferred to the intermolecular bond. 
Bond energies within both quantum models are not constant prior to a 
shock. Since energies are dependent upon mode population, total energy 
within a bond may oscillate. Energy may be transfered from bond to 


bond. The total energy of all bonds pnor to a shock should be constant. 
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(8) Deviations from the expected results of [Wyatt, Marston] in a harmonic 
model were caused primarily by the Weeeiee of states which resulted in 
normal oscillator motion of the atoms and center of mass motion to be 
inconsistent. Careful examination of the mode frequencies yielded 
information concerning the proper number of states one must have in 
order to describe atom motions properly. Lack of a substantial number of 
states also restricted the amount of energy transferred into the lattice. 

Upon re-examination of the problem, one sees that the total quantum 

Hamiltonian was written as follows: 

Hy = H° +V, +V(t) 
The anharmonic substituent (V4) of the Hamiltonian was formulated using a coupled 
normal mode operator: 

Ve= vm)? 264 + 60504) 

Instead of using a Hamiltonian which had all four modes coupled, we should have 
formulated the problem using a harmonic uncoupled Hamiltonian for modes one and 
two, and a coupled anharmonic Hamiltonian for modes three and four. 

Hy= Hy +H + H34+ Vaq+ V(t) + Volt) + V34(t) 

where 

Hry= total Hamiltonian 

H; = Harmonic Hamiltonian for mode 1 

H}4= Harmonic Hamiltonian of coupled modes three and four 

ie 4= Anharmonic Hamiltonian of coupled modes three and four 

V. (t) = Potential for mode 1 


V34(t) = Potenual for coupled modes three and four 








This would have allowed for an increased number of low frequency mode state (modes 
one and two); thus, placing the energy levels of all the modes in proper ratios. 

In conclusion, an entirely coupled Hamiltonian severely restricted proper 
formulation of this problem; however, reformulation using uncoupled modes should 
prove beneficial. Since desired results were not obtained using the harmonic model, 
deviations are most likely present in our anharmonic model. Once the Hamiltonian is 
uncoupled, an extremely large number of states may be added, and it will afford 


extremely large additions of energy to the crystal lattice. 





Table | 


Equilibrium values for position, atom separations, mass, 
and force constants 


qe = -6.9316 a, qo = +6.9316 a, 
qi = -3.4658 a, q3 = +1.7329 a, 
q> =-1.7329 a, q4 = +3.4658 a, 
pean 329 a, R, = +3.4658 a, 
my) = m3 = 34629 amu Ke = 0.00468 Hartree/Bohr2 
m2 = m4 = 34629 amu K = 0.14985 Hartree/Bohr? 


tb = 17314.5 amu 


X = -3.0 x 104 Hartree/Bohr? 





Table 2 


Frequencies and energy spacings for the normal modes. 


Normal Mode Frequencies (a.u.) 


Q 1 = 0.000256 
QO» = 0.000445 
€23 = 0.002927 


Q4 = 0.002949 


Energy Spacing (hQ)(eV) 


6.966 x 10-3 
12.108 x 10-3 
79.644 x 10-3 


80.243 x 10-3 
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Table 3 


Eigenvalues for the anharmonic and harmonic quantum 16 state system. 


State (n4 n3 n2 nq) Anharmonic Energy (a.u.) | Harmonic Energy (a.u.) 


OOOO 0.003284 0.003290 
0001 0.003540 0.003546 
0010 0.003722 0.003735 
OO11 0.003989 0039072 
0100 0.006163 0.006217 
1000 0.006235 0.006239 
0101 0.006435 0.006473 
1001 0.006594 0.006496 
0110 0.006634 0.006662 
1010 0.006684 0.006685 
ONG 0.006896 ORO eae 
1011 0.006942 0.006942 
1100 0.009165 0.009166 
1101 0.009421 0.009423 
IG) 0.009610 OO 
eed | 0.009859 0.009869 
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Figure 1 
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Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 
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Figure Captions 


Representation of the four atom cluster , with q< and as. 

Q< = initial cartisian coordinate position of entering qrcup atom 

q) = initial cartesian coordinate position of atom one 

q2 = initial cartesian coordinate position of atom two 

q3 = initial cartesian coordinate position of atom three 

q4 = initial cartesian coordinate position of atom four 

q> = initial cartesian coordinate position of the absorbing atom 

Ks = force constant utilized between diatom pairs 

K. = force constant utilized for entering, absorbing, and coupling of 
diatoms 

Ro = cartesian distance utilized for entering, leaving, and coupling of 
diatoms 


Mo = Cartesian distance utilized for coupling of diatom pairs 


Pictoral representation of the four normal modes. 

Normal mode one represents center of mass motion of the four atoms. 
Normal mode two represents paired diatom stretching about the central 
bond between atoms two and three. 

Normal mode three repesents an attractive stretching between atoms one 
and two, and a repulsive interaction between atoms three and four. 
Normal mode four represents a repulsive interaction between atoms one 
and two, and an attractive interaction between atoms two and three. 


Classical harmonic position of atoms qe, q1, 42, q3, q4, and qs. 
Positions are measured in atomic units, periods correspond to one 
harmonic oscillation of a diatom pair. Atom positions are measured 
from the center of mass for the chain. 


Classical anharmonic position of atoms de, q1, 42, 43, q4, and qs. 
Positions are measured in atomic units, periods correspond to one 
harmonic oscillation of a diatom pair. Atom positions are measured 
from the center of mass for the chain. 


Position of the ballistic particle. Atom position is measured in atomic 
units. Position is measured from the center of mass for the chain. 


Classical harmonic bond energy of the intramolecular bond (R12 ) 
between atoms q) and q2. Energy is measured in harmonic quanta. 


Classical harmonic bond energy of the intramolecular bond ( R34 ) 
between atoms q3 and q4. Energy is measured in harmonic quanta. 


Classical harmonic bond energy of the intermolecular bond ( R23 ) 
between atoms q? and q3. Energy 1s measured in harmonic quanta. 


Figure 9 
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Figure 12 


Figure 13 


Figure 14 


Figure 15 


Figure 16 


Figure 17 


Figure 18 


Figure 19 


Figure 20 


Figure 21 
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Classical anharmonic bond energy of the intramolecular bond ( R12 ) 
between atoms qj and q2. Energy is measured in harmonic quanta. 


Classical anharmonic bond energy of the intramolecular bond ( R34 ) 
between atoms q3 and q4. Energy is measured in harmonic quanta. 


Classical anharmonic bond energy of the intermolecular bond ( R23 ) 
between atoms q2 and q3. Energy is measured in harmonic quanta. 


Representation of the occupation of the two lowest energy states in the 
16 state quantum harmonic model. The probablities (nq n3 nz nq) 
correspond to the quantum number of a particular mode. 


Representation of the occupation of energy states (0010) and (0011) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (ng n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (0100) and (1000) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (ng n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (0101) and (1001) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (nq n3 n2 n}) of a particular mode. 


Representation of the occupation of energy states (0110) and (1010) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2 1) of a particular mode. 


Representation of the occupation of energy states (0111) and (1011) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (nq n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (1100) and (1101) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (ng n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (1110) and (1111) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (nq n3 n2 n}) of a particular mode. 


Representation of the occupation of the two lowest energy states (0000) 
and (0001) in the 16 state quantum anharmonic model. The numbers 
correspond to quantum numbers (n4 n3 n2 nq) of a particular mode. 


Representation of the occupation of energy states (0010) and (QO11) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2 nz) of a particular mode. 
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Figure 23 
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Figure 25 


Figure 26 


Figure 27 


Figure 28 


Figure 29 


Figure 30 
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Figure 33 


Figure 34 
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Representation of the occupation of energy states (0100) and (1000) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (0101) and (1001) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (0110) and (1010) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n?2 n1) of a particular mode. 


Representation of the occupation of energy states (0111) and (1011) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (nq n3 n?2 nq) of a particular mode. 


Representation of the occupation of energy states (1100) and (1101) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2 nj) of a particular mode. 


Representation of the occupation of energy states (1110) and (1111) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2 n}) of a particular mode. 


Representation of the ground state occupation of (OQOQOQ) in the 100 state 
quantum anharmonic model. The numbers correspond to quantum 
numbers (nq n3 n2 nj) of a particular mode. 


Representation of the occupation of the first excited state (0001) in the 
100 state quantum anharmonic model. The numbers correspond to 
quantum numbers (nq n3n72n)) of a particular mode. 


Representation of the occupation of excited energy state (0010) in the 
100 state quantum anharmonic model. The numbers correspond to 
quantum numbers (ng n3n92 n}) of a particular mode. 


Representation of the occupation of excited energy state (1134) in the 
100 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n4 n3 n2n}) of a particular mode. 


Representation of the occupation of excited energy state (1143) in the 
100 state quantum anharmonic model. The numbers correspond to 
quantum numbers (nq n3 n2 nj) of a particular mode. 


Representation of the occupation of excited energy state (1144) in the 
100 state quantum anharmonic model. The numbers correspond to 
quantum numbers (nq n3n2 nj) of a particular mode. 


Comparison of the positions of atoms one and two (qj,q2) in the 16 
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Figure 36 


Figure 37 


Figure 38 


Figure 39 


Figure 40 


Figure 41 


Figure 42 


Figure 43 


state and 100 state quantum harmonic models. Postions are measured in 
atomic units (ag) and one period corresponds to one oscillation of an 
atom. - 


Comparison of the positions of atoms three and four (q3,qq4) in the 16 
state and 100 state quantum harmonic models. Postions are measured in 
atomic units (a9) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms one and two (qq,q2) in the 64 
state and 100 state quantum harmonic models. Postions are measured in 
atomic units (a9) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms three and four (q3,q4) in the 64 
state and 100 state quantum harmonic models. Postions are measured in 
atomic units (ag) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms one and two (q},q2) in the 16 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (ag) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms three and four (q3,q4) in the 16 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (ag) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms one and two (q},q2) in the 64 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (ag) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms three and four (q3,qq4) in the 64 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (a9) and one period corresponds to one oscillation of an 
atom. 


Comparison of the positions of atoms one and two (qj,q2) in the 100 
state quantum harmonic model and the classical harmonic model. 
Postions are measured in atomic units (a 9) and one period corresponds 
to one oscillation of an atom. 


Comparison of the positions of atoms three and four (q3,q4) in the 100 
state quantum harmonic model and the classical harmonic model. 
Postions are measured in atomic units (a9) and one period corresponds 
to one oscillation of an atom. 
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Figure 44 


Figure 45 


Figure 46 


Figure 47 


Figure 48 


Figure 49 


Figure 50 
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Figure 52 
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Figure 55 


Figure 56 
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Comparison of the positions of atoms one and two (q},q2) in the 100 
state quantum anharmonic model and the classical anharmonic model. 
Postions are measured in atomic units (a9) and one period corresponds 
to one oscillation of an atom. 


Comparison of the positions of atoms three and four (q3,q4) in the 100 
state quantum anharmonic model and the classical anharmonic model. 
Postions are measured in atomic units (a9) and one period corresponds 
to one oscillation of an atom. 


Varations in atomic position (a9) of atom one (q1) in the quantum 100 
state anharmonic and harmonic models. 


Varations in atomic position (a9) of atom two (q2) in the quantum 100 
state anharmonic and harmonic models. 


Varations in atomic position (a9) of atom three (q3) in the quantum 100 
state anharmonic and harmonic models. 


Varations in atomic position (a) of atom four (q4) in the quantum 100 
state anharmonic and harmonic models. 


Positions of atoms q and q? in the 16 state harmonic model in which 
the potential created by qs is zero at all ttmes. Positions are measured in 
atomic units. : 


Positions of atoms q3 and qq in the 16 state harmonic model in which 
the potential created by q, is zero at all times. Positions are measured in 
atomic units. 


Positions of atoms q; and q? in the 16 state harmonic model ; the 
potential created by qs is applied as a step jump at t = 10 periods. 
Positions are in atomic units. 


Positions of atoms q3 and qq 1n the 16 state harmonic model ; the 
potential created by qs is applied as a step jump at t = 10 period. 
Positions are measured in atomic units. 


Positions of atoms q; and q? in the 16 state harmonic model ; the 
potential created by qs is applied as a ramp increase between t = 10 and 
12.5 periods. Positions are measured in atomic units. 


Positions of atoms q3 and qq in the 16 state harmonic model ; the 
potential created by q; is applied as a ramp increase between t = 10 and 
12.5 periods. Positions are measured in atomic units. 


Positions of q, and q2 in the 16 state harmonic model ; the potential 
created by qs is the actual potential of qs shifted 3 pernods earlier. 
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Figure 63 


Figure 64 


Figure 65 


Figure 66 


Figure 67 


Figure 68 
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Positions are measured in atomic units. 


Positions of q3 and qq in the 16 state harmonic model ; the potential 
created by qs is the actual potential of q, shifted 3 periods earlier. 
Positions are measured in atomic units. 


Positions of qy and q? in the 16 state harmonic model ; the potential 
created by qs is the actual potential of qs shifted 3 periods later. 
Positions are measured in atomic units. 


Positions of q3 and qq in the 16 state harmonic model ; the potential 
created by qs is the actual potential of qs shifted 3 periods later. 
Positions are measured in atomic units. 


Positions of atoms qj and q? in the 100 state quantum anharmonic 
models. The solid line plots are atom positions with gamma = -3.0 x 10 
-4 Hartree/ Bohr 3. The dashed line plots indicate atom positions with 
ten times the anharmonicity. 


Positions of atoms q3 and qq in the 100 state quantum anharmonic 
models. The solid line plots are atom positions with gamma = -3.0 x 10 
-4 Hartree/ Bohr 3. The dashed line plots indicate atom positions with 
ten times the anharmonicity. 


Distribution of energy into mode one of the 100 state quantum harmonic 
model. Energy is measured in quanta. 


Distribution of energy into mode two of the 100 state quantum harmonic 
model. Energy is measured in quanta. 


Distribution of energy into mode three of the 100 state quantum 
harmonic model. Energy is measured in quanta. 


Distribution of energy into mode four of the 100 state quantum 
harmonic model. Energy is measured in quanta. 


Distribution of energy into mode one of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 


Distribution of energy into mode two of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 


Distribution of energy into mode three of the 100 state quantum 
anharmonic model. Energy is measured 1n quanta. 


Distribution of energy into mode four of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 


Potential and kinetic energy of bond ( R12 ) between atoms q, and q? in 
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Figure 76 


Figure 77/ 


Figure 78 


Figure 79 


the 100 state quantum harmonic model. Energy is measured in quanta. 


Potential and kinetic energy of bond ( R34 ) between atoms q3 and qg in 
the 100 state quantum harmonic model. Energy is measured in quanta. 


Total energy of bonds (R12 ) and ( R34 ) in the 100 state quantum 
harmonic model. Energy is measured in quanta. 


Potential and kinetic energy of bond ( R23 ) between atoms q? and q3 in 
the 100 state quantum harmonic model. Energy is measure in quanta. 


Total energy of bond (R23 ) between atoms q? and q3 in the 100 state 
quantum harmonic model. Energy is measured in quanta. 


Potential and kinetic energy in bond ( R12 ) between atoms q, and q? in 
the 100 state quantum anharmonic model. Energy is measured in 
quanta. 


Potential and kinetic energy in bond ( R34 ) between atoms q3 and q4 in 
the 100 state quantum anharmonic model. Energy is measured in 
quanta. 


Total energy in bonds (R12 ) and ( R34 ) in the 100 state quantum 
anharmonic model. Energy is measured in quanta. 


Potential and kinetic energy of bond ( R23 ) between atoms q? and q3 in 
the 100 state quantum anaharmonic model. Energy is measured in 
quanta. 


Total energy in bond (R23 ) between atoms q? and q3 in the 100 state 
quantum anharmonic model. Energy is measured in quanta. 
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